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NOMENCLATURE 


a 



H 

k 


n 
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(r , 9 , z) 

(r*,0*,z*) 

R i 

R(z) 

(U,V,W) 
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w i 
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speed of sound 

group velocity of wave propagation In axial direction; -du/dk 

phase velocity of wave propagation In axial direction; -w/k 

total head function 

axial wave number 

azimuthal wave number 

pressure of basic flow field 

pressure of disturbance flow field 

amplitude function of pressure disturbance 

ratio of swirl velocity to axial velocitv at inlet of divergent 
duct 

cylindrical coordinates at inlet of duct; r and z are normalized 
by R t 

local cylindrical coordinates; r* - r/R(z s ), 0* = 0, z* - z - z s 
radius at duct inlet 

geometric shape of duct normalized by Rj 

velocities of basic flow field; U in radial direction, W in 
axial direction 

velocities of disturbance flow field; u in radial direction, 
w in axial direction 

amplitude functions of velocity disturbances; 4/ in radial direc- 
tion, W in axial direction 

uniform axial velocity at inlet of divergent duct 

stream function of basic flow field 

circulation function 

constant angular velocity 
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frequency of disturbance 


p water density 

Subscripts 

c critical state 

i function evaluated at lr.let of divergent duct, z ■ 0 

o function evaluated at r ■ 0 

s axial station 
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AN EXAMINATION OF A GROUP-VELOCITY CRITERION FOR THE BREAKDOWN 


OF AN IDEALIZED VORTEX FLOW 
Chon- Yin Tsai* 

Ames Research Center 
and 

Sheila E. Widnall** 
Massachusetts Institute of Technology 


SUMMARY 


The phenomenon of vortex breakdown is believed to be associated with a 
finite-amplitude wave that has become trapped at the critical or breakdown 
location. The conditions at which the propagating waves become trapped at a 
certain axial location is examined here by use of a group-velocity criterion 
implied by Landahl's general theory of wave trapping. Contrary to previous 
theories, the present explanation may cover both axisvmmetric and asymmetric 
breakdown. In order to investigate the relationship between wave trapping and 
vortex breakdown, a numerical study is made here of an ideal vortex having 
constant vorticity and uniform axial velocity at the inlet of a slowly diverg- 
ing duct. The solution presented by Batchelor for the flow field of a vortex 
flow in a slowly varying duct is well known. The linear wave propagation 
analysis is then applied to this base flow at several axial stations for sev- 
eral values of the ratio of swirl velocity to axial velocity at the inlet of 
the divergent duct, assuming a locally parallel flow. The dispersion rela- 
tions and hence the group velocities of both the symmetric (n = 0) and asym- 
metric modes (n = ±1) are investigated. The existence of a critical state in 
the flow (at which the group velocity vanishes), and its relationship to the 
stagnation point on the axis of the duct and to the occurrence of an irregular 
singularity in the equations governing wave propagation in the flow field are 
discussed. 


INTRODUCTION 


The study of the stability of vortex flows is important in many fields, 
such as aeronautics, combustion, and geophysics. One of the instability 
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mechanisms of vortex flows is the so-called vortex breakdown or vortex burst- 
ing. That is, the breakdown of vortex motion is identified as occurring when 
the swirl and axial flow undergo an abrupt and drastic structure change along 
the axi3 of the vortex. The breakdown or burst is often characterized by the 
tormation of an internal stagnation point on the vortex axis (ref. 1), fol- 
lowed by a reversed flow in a region of limited axial extent. The first 
documented vortex breakdown phenomenon was reported by Pockham and Atkinson 
(unpublished report, 1957) in an investigation of the aerodynamics of delta 
wings. Vortex breakdown also has been observed in the swirling flows through 
nozzles and diffusers. Practical applications of the vortex breakdown 
phenomenon occur in two ways. In some cases, one wishes to avoid or delay the 
occurrence of vortex breakdown, as in the flow over a delta wing, to control 
stall. In other cases, one wishes to Induce vortex breakdown to dissipate 
strong velocity gradients in the trailing vortices behind large aircraft during 
takeoff and landing to reduce the hazard to following aircraft or to reduce 
blade/vortex interaction noise of helicopter rotors. Solutions to ..hese prob- 
lems require an improved understanding of the vortex breakdown phenomenon. It 
is important to know which vortices will suffer breakdown and what are the 
malor features of the recirculation zone and of the downstream flow field. 

In order to Isolate the phenomenon under more controlled circumstances, 
experiments on vortices are usually performed in ducts. Harvey (ref. 2) used 
a cylindrical vortex formed in a long duct of constant cross section to study 
vortex breakdown. By varying the amount of swirl imparted to the fluid before 
it entered the duct, he found that the breakdown was the intermediate stage 
between two basic types of rotating flows. He concluded that breakdown is a 
characteristic of "critical" phenomenon and that it is not a hydrodynamic 
instability. Sarpkaya (ref. 3) obtained breakdown in a divergent duct with 
less swirl than that required in a duct of constant cross section. He observed 
three types of vortex breakdown: double helix (mild), spiral (followed by 

turbulent mixing), and axisymmetric (followed by a thicker vortex core, then a 
spiral breakdown, and, finally, by turbulent mixing). Faler and Lebovich 
(ref. 4) conducted experiments in an apparatus modeled after Sarpkaya's and, 
from their flow visualization studies, found three additional states of vortex 
disruption. Their experimental results show that the type and location of 
breakdown were dependent on the Reynolds number and circulation distribution 
of the flow. The first detailed velocity measurements of the vortex flow 
field containing breakdown in a slightly diverging duct were obtained by using 
a laser Doppler anemometry (refs. 4-6) . Although theoretical attempts to 
explain the vortex breakdown phenomenon are numerous, none have led to a sat- 
isfactory understanding of the observed flow. The various theories for vortex 
breakdown have been reviewed and classified into three categories by Hall 
(ref. 7) and Lebovich (ref. 1): (a) vortex breakdown is proposed to be the 

consequence of hvdrodvnamic instability, (b) vortex breakdown is analogous to 
the separation of a boundary laver in an adverse pressure gradient, and 
(c) vortex breakdown is the result of some type of critical condition existing 
in the flow. Theories using hydrodynamic instability as the sole mechanism 
for vortex breakdown are no longer considered satisfactory (e.g., Ludwieg's 
theory (ref. 8) is mathematically inconsistent if applied to duct flow (see 
review by Lebovich (ref. 1)). The analogy between vortex breakdown and the 
separation of a boundary layer does not suggest any physical mechanism for the 
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phenomenon and is incomplete because of the restriction to axisvmmetric flow. 
The idea of a critical condition existing in a vortex flow was Introduced bv 
Squire (ref. 9). Based on phase speeds. Squire proposed that if infinitesi- 
mal amplitude axisvmmetric waves have zero phase speed, then the disturbances 
originating downstream may propagate upstream and cause breakdown. Although 
Squire's explanation of vortex breakdown has some weaknesses, his idea of the 
possible appearance of infinitesimal standing waves in the flow has led to the 
development of various wave theories for vortex breakdown. Benjamin (ref. 10) 
considered vortex breakdown as precisely analogous to the hydraulic jump in 
open channel flow. That is, the frictionless swirling flows generally occur 
in conjugate pairs wherein both can form parts of the same overall system. 

For a given distribution of total head and circulation over the stream sur- 
faces, one possible state of flow is defined as "subcrlt ical" — an Infinitesi- 
mal axlsymmelric standing wave can occur. Its conjugate state is defined as 
"supercritical" — no standing wave is possible. The analytical condition for 
the existence of infinitely long, axisynunetric standing waves, Cp ■ 0, has a 
fundamental role in Benjamin's theory. 

Randall and Lebovich (ref. 11) developed a theory of long, weakly nonlin- 
ear waves propagating along critical flows in tubes of variable cross section. 
Although the structural features of the vortex breakdown predicted bv a non- 
linear trapped wave model are consistent with manv of the experimentally 
observed features, there are still deficiencies in the theory — reversed 
swirls predicted in the bubble interior are physically impossible and the 
spiral form of breakdown is excluded. Numerical experiments with the Navier- 
Stokes equations restricted to axial symmetry and steady flow throughout the 
entire flow field, including the breakdown region, have been studied by sev- 
eral investigators. Although there are some qualitative points of agreement 
between physical and numerical experiments, there are a number of ways in 
which the numerical and physical experiment are quite different (see ref. 1). 
The most common deficiencies in previous theories are that the spiral form of 
breakdown is excluded and the recirculation zone itself is assumed to be 
axisymmotric. Therefore, a different approach to the study of vortex break- 
down is required. 

Using Landahl's (ref. 12) general theory of wave focussing as the frame- 
work for analysis, Bilanin (ref. 13) suggested and explored the possibility 
that vortex breakdown is associated with a trapped finite-amplitude wave. He 
defined the vortex flow as "supercritical" or "subcritical" when the group 
velocity of infinitesimal waves is in the "downstream" or in "both" directions 
along the axis of the vortex. Thus, the same flow may be subcritical to axi- 
svmmetric waves and supercritical to spiral wave modes. He then obtained the 
critical state (defined as occurring when C_ - 0) for the propagation of 
various wave types along a line vortex undergoing solid body rotation with 
uniform axial velocity. By ranking the modes according to their ability to 
signal upstream, he postulated which type of breakdown (axisvmmetric vs. 
spiral) would first occur. 

In the present study, the linear wave trapping process is explored for an 
ideal model of vortex flow — for example, a vortex having constant vorticitv 
n and constant axial velocity Wj at the inlec of a slowlv diverging duct . 
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This flow model has the simplest entry conditions and provides a hypothetical 
configuration that has a minimum number of variable parameters. Although this 
flow is not a good model of vortex breakdown because it has negative swirls 
downstream of the stagnation point, the wave characteristics of these profiler, 
may be of some Interest. 

In the following, we first discuss the group-velocity criterion for vortex 
breakdown as derived from Landahl's (ref. 12) general wave trapping theory. 
Then, the lineai wave propagation analysis is formulated and applied to the 
idealized model of vortex flow in a slightly diverging duct obtained in ana- 
lytic form by means of the invlscid quasi-cyl lndrlcal approximation. The dis- 
persion relations and hence the group velocities are then calculated 
numerically. 

Dr. Vernon J. Kossow suggested the application of wave trapping theory to 
the theoretical model of uniform rotation and axial velocity. The authors 
thank Dr. Rossow for valuable discussions and encouragement. 


WAVE TRAPPING THEORY AM) VORTEX BREAKDOWN 


Landahl (ref. 12) developed a general theory for the wave mechanics of 
breakdown to determine under what conditions the steady or unsteady laminar 
flow will break down into high-frequency oscillations. The analysis of a 
small-scale secondary wave riding on a large-scale inhomogeneity reveals that 
the breakdown me hsnism occurs when the group velocity of a self-excited sec- 
ondary wave is near the phase velocity of the primary wave. To Illustrate, a 
simple example of wave trapping given by Landahl is included here. Consider 
the supersonic one-dimensional gas flow in a Laval nozzle. An oscillating 
disturbance produces two sets of acoustic waves: an advancing wave with a 
propagation velocity W + a (where W and a are the local flow and sound 
velocity, respectively) and a receding wave with a propagation velocity W - a. 
Linear theory predicts the amplitude of the receding wave as (W - a) -1 ' 2 . 
Therefore, if sonic flow exists at the throat, the amplitude of the waves 
would become large and could be predicted only by a nonlinear analysis. The 
singularity in the linear theory is a result of a concentration of wave energy 
produced by consecutive wave fronts catching up with each other. The wave 
trapping process is a possibility for any inhomogeneous continuum system which 
can support propagating waves that are small relative to the scale of inhomo- 
geneities. The large-scale inhomogeneity could be produced by a primary wave, 
by viscous dissipation, or by changes in the boundaries of the flow field. In 
the present study, the inhomogeneity is assumed to be induced bv the slowly 
diverging wall of the duct. Landahl's breakdown criterion for steady vortex 
flow In a duct requires that the group velocitv, C R , of infinitesimal waves be 
zero at some noint. The vortex flow mav then be classified as ’’sunercrit leal" 
to a wave type when the group velocity of that wave propagates in the down- 
stream direction (Cg >0), and "suberit ical" if the wave can propagate in 
either direction. The dispersion relations and hence the group velocities are 
obtained by solving the eigenvalue problem at each axial station as if the 
flow were locally parallel. The location where the group velocity of wave 
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disturbance Is zero Is proposed as the place where wave trapping and therefore 
the breakdown of vortex flow would occur. Typically, only one type of wave 
mode (axlsymmctric or spiral) would be critical at any station in the flow, 


FORMULATION OF THE PROBLEM 


It was proposed in the previous section that the phenomenon of vortex 
breakdown Is associated with the trapping of f ini te-amnlitude waves propagat- 
ing along the vortex. To investigate the relationship between wave structure, 
wave trapping, and vortex breakdown, we consider a slowlv varying basic flow 
field perturbed by axisvrametric and asymmetric disturbance waves. The problem 
is formulated in the following two sections. 


Basic Flow Field 


In the present study, the basic flow field is taken as a swirling flow 
inside a duct with a small divergent angle. Using cylindtical polar coordi- 
nates (r,6,z) with the corresponding velocity components (U,V,W) and intro- 
ducing the stream function 'Hr.z) such that U • (-1/r) (3^/3z) and 
W « (1/r) (d<P/dv ) , the inviscid, incompressible equations of motion can be 
reduced to (ref. 14) 


where 


af* af* _ I „ 2 dH _ r — 

3z 2 3r 2 r 3r r d<|/ d>J 


H(H 


1 (U 2 + V 2 + W 2 ) + — (Total head function) 

2 p 

r('I’) ■ rV (Circulation) 


( 1 ) 


The pressure P is calculated from the radial equation of motion. 


For the relatively simple case considered here, the fluid at the inlet of 
the divergent duct has uniform axial velocity Wj and rotates as a rigid body 
with angular velocity fl; the entrance conditions are therefore given bv 


and 


r - flri 2 « ♦ 

i w t 


dH __i iH li* 

d* ‘ Wi ' r 2 j 2 ' Wj 2* 3 2 
l z z 
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where Che subscript 1 of r^ denotes the radial coordinate at the t'let of 
the divergent Juct. The dependence of H and T on must hold throughout 
the flow field. The governing equation (1) then takes the form 


a 2 * . a 2 * 1 a* 2n 2 2 4fl 2 . , r 2 a 2 ^ 

IP + IP - r li * ^ 1 - * * H 1 2? IP 


( 2 ) 


For the flow in a slowlv diverging duct, the quasi-cvlindrical approximation 
is applicable because the vaiiations in axial direction are much smaller than 
variations in the radial diruction. Therefore, equation (2), in dimensionless 
form with length scale Rj (radius at the duct inlet) and velocity scale Wj, 
becomes 


iii _ i ut 

3r 2 " r 3r 


2q 2 r : - 4q 2 <|< 


The boundary conditions are 


MO) - 0 


MR) 


1 

2 


where R is the dimensionless local radius of the linearly divergent duct and 


r a 

q " w 1 

The stream function for the flow field is written in analytic form as (ref. 7) 

* - | r 2 + \ F(R)rJ!(2qr) (3) 


with 

F(R) “ RJ i (2qR) 

where Jj is the Bessel function of the first kind, of order 1. The quasi- 
cylindrical solution (3) fails whenever J i ( 2qR) » 0. This failure of the 
quasi-cylindrical approximation has often been taken as coincident with the 
occurrence of vortex breakdown. It will be seen from the results obtained in 
the next section that C (for mode n - 1) tends to zero at an axial location 
much earlier than the location where the quasi-cylindrical approximation fails 

From the stream function, the axial and swirl velocities can be calcu- 
lated and written as: 


V = 2q ^ - q [r + F(R)Jj(2qr)] 


(4) 
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( 5 ) 


1 ail* 

w - 7-57 - 1 ♦ qF(R)J 0 (2qr) 

where 0 i r £ R and J Q is the 1 es 3 el function of the first kind, of order 0. 
The streamline patterns for various q in a duct of shape function R(z) are 
given in figures 1 through 5. The z term (normalized by the inlet radius 
Rj) is the downstream axial coordinate starting at the duct inlet. It can be 
seen that the reversed swirl velocities occur Inside the streamline ^ ■ 0 
starting from the stagnation point. Since T (V' ) is undefined for < 0, the 
solution is meaningless beyond the stagnation point. 


Linear Wave Propagation Analysis of a Swirling Flow 

The wave propagation analysis of a rotating flow was first studied by 
Thomson (ref. 15). He investigated the vibration of an isolated column vortex 
with solid-body rotation surrounded by an irrotational revolving fxow and 
showed that it is stable to all infinitesimal disturbances. Recently, analy- 
ses have been made of the stability of flow fields wherein a sheared flow is 
superimposed on a rotating flow (e.g., refs 16-20). As an application of 
these analyses, the propagation of an inviscid linear wave along a swirling 
flow (eqs. (4) and (5)), bounded by a rigid cylindrical duct 0 i r* i 1 is 
investigated, where r* «• r/R. The coordinates (r*,0*,z*) for the wave analy- 
sis is at ee :h axial station along the duct, assuming that the base flow is 
locally parallel. This basic flow field is subjected to small perturbations 
which permits the equations of motion to be linearized. The small perturba- 
tions are assumed to be of the form: 

(u, \>, w, p/p } = (3f(r*) , y(r*) t , ^(r*)} exo[i (u>t + n9* + kz*) ] 

so that the group velocity is given by 

C ■ -dw/dk 

and the phase velocity is given by 


C p - -u./k 

where u, \>, w, and p are the perturbation quantities of radial velocity, 
tangential velocity, axial velocity, and pressure, respectively. (We drop 
superscript * in the following analysis for convenience.) The resulting single 
equation for W (r) is 


d ? W 

dr 2 



(— 

dS 

+ 1 

da B 

- 3&X 

/ dr 

V>s 

dr 

Y 

dr y 2 s 

Y 2 / 


4/ 


( 6 ) 


where y * u + nV/r + kW, k and .1 are axial and azimuthal wave numbers, 
to ■ w r + iii)^ is the complex frequency, ana 
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The boundary condition at r - 0 is 


“*( 0 ) 


0 if n i ±1 

Finite if n ■ ±1 


The boundary condition at the wall, r ■ 1, is 

♦( 1 ) - 0 


(7) 


( 8 ) 


Equation (6) along with the boundary conditions (7) and (8) constitute 
an eigenvalue problem. Instability will result if the imaginary part of the 
eigenvalue ui is negative for any given pair of wave numbers (n,k). 

For the axisymmetric standing wave (i.e., n ■ 0 and w ■ 0) , which is the 
case considered by Benjamin (ref. 10), equation (6) is written as 


d 2 % Id ty _ f 1 1 / d 2 W 1 dW\ _2V_ /dV V\ 2 

dr 2 r dr "[ r 2 w \dr 2 " r dr / “ rW \dr r/ 


* - 0 


(9) 


If we -ibstitute the inviscid quasi-cylindrically approximate solution (4) 
and ^j) for the basic flow field, equation (9) can be reduced to the form: 


0 + T 3? * (? - 4 ’ 2 + “ 2 )* * 0 <10 > 

which is independent of the geometry of duct R. 

The solution to equation (10), which is well behaved at r ■ 0, is 

'fc'(r) - AJ ! [ (4q 2 - k 2 ) l/2 r] 

The boundary condition at r * 1 gives 

k 2 - 4q 2 - b n 2 

•./here b n ■ 0, 3.85, 7.02, . . . are zeros of the Bessel function Ji. 
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According to Benjamin's (ref. 10) conjugate flow theory, a given state of 
flow is supercritical if all eigenvalues k‘ are negative, so that no stand- 
ing wave ir possible. The disturbances will then diminish in the vpstream 
direction as exp(-kz). The state of the flow is considered subci itical if at 
least the first of k is positive. Therefore, the critical state for the 
special case considered here occurs when q c - b n /2 * 1.92. That is, if the 
ratio q of the vortex flow with solid-body rotation and uniform axial veloc- 
ity at the inlet of a divergence duct is less than q c , the flow is supercriti- 
cal. Since equation (10) is independent of duct, geometry, the state of the 
flow in the whole duct is always supercritical according to the Benjamin's 
theory and it will not be driven to a subcritical state by the divergence of 
the wall. This, however, is not the case if the present definition of criti- 
cal state by the group-velocity criterion is used. 

Except for the special case of uniform vortlcity and axial flow (inlet 
condition), the dispersion relation for a general profile must be obtained 
numerically. Hultgren (ref. 21) presented a numerical method which uses 
Moulton's method for the stability calculation of rotating gas flows. By use 
of this method, it was possible to obtain the eigenvalues for the velocity 
profiles (4) and (5) at each axial station for various q values and deter- 
mine the group velocity, assuming a locally parallel flow. 

Equation (6) has a regular singularity at r - 0. A Frobenius power 
series solution is used to start the numerical Integration at the axis of the 
duct. The solution is then advanced toward the boundary at r • 1. Another 
regular singularity may occur at some value of r in the range of 0 < r s 1 
tor which y(r) ■ 0 (y(0) ■ 0 gives an irregular singularity at r ■ 0 and 
must be handled differently). For y(r) “0 in the range 0 < r i 1, the 
integration path must be deformed in agreement with Lin's (ref. 22) criterion. 
The flows of interest (i.e., on which the group velocity C g has a minimum 
value) in the present problem turn out to be outside the range y(r) ■ 0 for 
0 < r & 1. 

The dispersion relations and hence the group velocities for axisymmetric 
(n ■ 0) and asymmetric (n * ±1) waves have been obtained numerically for the 
mean profiles at each axial station considered as a parallel flow and are pre- 
sented in the next section. 


RESULTS AND DISCUSSION 


In this section, the numerical results for the characteristics of wave 
propagation along the duct are presented. The velocity distribution in 
duct is taken as the solution found by the inviscid quasi-cylindrical approxi- 
mation. The analysis is applied at various axial stations along the divergent 
due for q ■ 0.8, 1.0, 1.429, 1.667, and 2.0 (figs. 1-5). 

The group velocities for the axisymmetric mode (n * 0) are shown in the 
figures 6 through 1C. It can be seen that the group velocity for the flow 
region ahead of the stagnation point is directed downstream (i.e., C g > 0 for 
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all wave number k) and hence Che flow in supercritical ahead of Che stagnation 
point, except for q - 2.0. For the region behind the stagnation point, the 
group velocity is directed upstream and hence the flow is subcritical behind 
the stagnation point. Recall that the flow downstream of the stagnation point 
is physically meaningless in that T is undefined for negative ♦ . The wave 

propagation results are presented in this region of Interest. The slope of 

the amplitude function of radial perturbation velocity near r ■ 0 becomes 
larger as its axial location approaches the stagnation point. This occurs 
when y 0 ■ u + n(V/r)| r-0 ♦ kWl r . Q goes to zero at r ■ 0 so that the 
Frobenius power series solution near r ■ 0 (obtained from the eq. (6)) 
becomes invalid. 

Therefore, for the axisvmmetric mode (n ■ 0) of wave propagation, near 
the stagnation point of the basic flow field, the approach to the irregular 
singularity Y 0 * 0 causes the group velocity to go to zero. However, this 
observation is only qualitatively true since the parallel flow assumption Is 
not valid near the stagnation point. 

For the asymmetric mode (n ■ 1), the group velocities are presented in 

figures 11 through 15. Note that the group velocity generally decreases as 

the axial station moves downstream toward the stagnation point for supercriti- 
cal flows, except when q ■ 1.429 (fig. 13). Once again, C- may cross zero 
in the range 1.36 < R(z) < 1.43 and 1.2 < R(z) < 1.33 for q ■ 0.8 and 1.0, 
which is the region containing the stagnation point. If we rb ~k a typical 
amplitude function of the radial perturbation velocity (sho; i figs. 16 
and 17), we see that the slope of the amplitude functions which C_ is 
minimum, is very large near r - 0 as its axial station approaches the stagna- 
tion point. This situation is similar to the e> Asymmetric case discussed pre- 
viously. Figure 13 shows the group velocity at various axial stations for 
q - 1.429. It is seen that, at the inlet of the divergent duct, Cg mln is 
very close to zero, but the amplitude function (fig. 18) of the radial per- 
turbation velocity is well behaved. (Figure 19 shows the amplitude function 
at axial station behind the stagnation point.) As the station being studied 
moves downstream, Cg . does not cross zero, but increases as its axial sta- 
tion moves closer to the stagnation point. This result differs from the two 
cases for q ■ 0.8 and 1.0 discussed previously. For q ■ 1.667, the flow at 
the inlet of the duct is still supercritical to the axisymmetric waves 
(fig. 9), but is subcritical to the asymmetric (n ■ 1) waves (fig. 14). For 
q - 2.0 (which Is only slightly larger than Benjamin's (ref. 11) critical 
value of 1.92), figures 10 and 15 show that the flow is subcritical to both 
axisymmetric and asymmetric waves. However, figure 10 shows that this flow is 
subcritical to a mode having the simple structure of the amplitude function, 
but is still supercritical to a different mode having several nodes along the 
radial axis in its amplitude function. 

The group velocities of asymmetric mode (n ■ -1) at various axial sta- 
tions for various q values are shewn in figures 20 through 24. The flows 
are supercritical at the inlet of the duct except for q - 2.0. (Figure 24 
shows that C g at R(z) “ 1.0 crosses zero in the long wavelength limit.) 
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The numerical results for phase velocities are presented in figures 25 
through 38. The directions of group velocities and phase velocities and the 
states of the i lows are compared in table 1. For q £ 1.667 (which is less 
than q c ), the directions of group velocity and phase velocity are the same 
for modes n ■ 0 ai d l at the same axial locations. The exception is mode 
n ■ 1 for j 1.667. For mode n - -1 of all cases, and modes n • 0, il 
of case q ■ 2.0, the directions of the group velocity and phase velocity are 
not the same in a certain range of wave number k. 


CONCLUSIONS 


The linear wave propagation analysis has been applied to an idealized 
model of the flow in a slowly diverging duct having inlet conditions of solid- 
body rotation and uniform axial velocity. For a varying q (the ratio of 
swirl velocity to axial velocity at the inlet), the wave propagation of vari- 
ous modes at each axial station (treated as a parallel flow) has been calcu- 
lated. Although the development of the flow field in the divergent duct 
c. scribed by the lnviscid quasi-cylindrical approximation is unrealistic in 
the region near the downstream stagnation point, it may be of interest to 
identify the wave propagation characteristics of these different flow regions. 

At the present stage, we conclude the following for this idealized flow 
mode 1 : 

(a) If the critical state is defined as done by Benjamin (ref. 10), the 
flow that is supercritical cannot be driven to be subcritical simply by the 
divergent wall of the duct. However, if the critical state is defined by the 
group-velocity criterion Cg * 0, the flow that is supercritical can be driven 
to be subcritical due to the occurrence of the stagnation point by tbs diver- 
gent wall of the duct. 

'’b) As the stagnat* an point is approached, the group velocity C R goes 
to ze o for the ax i symmetric waves. However, the parallel flow assumption 
near the stagnation point is not valid. 

(c) The directions of phase velocity and group velocity are consistent 
only for modes n = 0 and 1 in the flow that is supercritical upstream of the 
stagnation point. 
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TABLE 1.- COMPARISON OF DIRECTIONS OF CROUP VELOCITY C g AND PHASE VELOCITY 
C p AND THE STATE OF THE FLOW CONDITION AHEAD OF THE STAGNATION POINT 


q 

Mode n • 0 

Mode n ■ 1 

Mode n ■ -1 

0.8 

Same 

Supercritical 

Sami’ 

Supercritical 

Not the same 
Supercritical 

1.0 

Same 

Supercritical 

Same 

Superct ltical 

Not the same 
Supercritical 

1.429 

Same 

Supercritical 

Same 

Supercritical 

Not the same 
Supercritical 

1.667 

Same 

Supercritical 

Not the same 
Subcritlcal 

Not the same 
Supercritical 

2.0 

Not the same 
Subcritlcal 

Not the same 
Subcritlcal 

Not the same 
Subcritlcal 
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a divergent duct for q = 1.0. 
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Figure 3.- Stream function in a divergent ^uct for q 
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Figure 4.- Stream function in a divergent duct for q 


Figure 6 


8 


6 


4 


c 


9 


2 


0 


-.2 
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- Group velocities of mode n = 0 at 
q - 0.8 at the inlet of a dive 
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Figure 8.- Group velocities of mode n * 0 at various axial stations for 
q = 1.429 at the inlet of a divergent duct. 


Figure 9 



- Group velocities of mode n ■ 0 at various axial stations for 
q = 1.667 at the inlet of a divergent duct. 
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Figure 12.- Group velocities of mode n - 1 at various axial stations for 
q - 1.0 at the inlet of a divergent duct. 



k 

Figure 13.- Group velocities of mode n - 1 at various axial stations for 
q * 1.429 at the inlet of a divergent duct. 
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Figure 14.- Group velocities of mode n • 1 at various axial stations for 
q - 1.667 at the inlet of a divergent duct. 
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Figure 15.- Group velocities of mode n = 1 at various axial stations for 
q - 2.0 at the inlet of a divergent duct. 




Figure 16.- Amplitude functions of radial perturbation velocities (mode n ■ 1) 
in a region ahead of the stagnation point for q - 0.8 at the inlet of a 
divergent duct. 



Figure 17.- Amplitude function of radial perturbation velocities (mode n ■ 1) 
in a region behind the stagnation point for q * 0.8 at the inlet of a 
divergent duct. 
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Figure 18.- Amplitude fu.i.;tlcns of radial perturbation velocities (mode n ■ 
in a region ahead of the stagnation point for q ■ 1.429 at the inlet of 
a divergent duct. 
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Figure 19.- Amplitude functions of radial perturbation velocities (mode n - 
in a region behind the stagnation point for q - 1.429 at the inlet of a 
divergent duct. 




Figure 20.- Group velocities 01 mode n - -1 at various axial stations for 
q ■ 0.8 at the inlet of a divergent duct. 



k 


Figure 21.- Group velocities of mode n ■ -1 at various axial stations for 
q » 1.0 at the inlet of a divergent duct. 
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Figure 22.- Croup velocities of mode n ■ -1 at various axial stations for 
q - 1.429 at the inlet of a divergent duct. 
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Figure 2J.- Group velocities of mode n - -1 at various axial stations for 
q - 1.667 at the inlet of a divergent duct. 
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Figure 24.- Group velocities of mode n * -1 at various axiaJ. stations for 
q ■ 2.0 at the inlet of a divergent duct. 
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Figure 25.- Phase velocities of mode n = 0 at various axial stations for 
q = 0.8 at the inlet of a divergent duct. 
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1 .;ure 27.- Phase velocities of mode n - -1 at the various axial -Hat ions 
for q ■ 0.8 at the inlet of a divergent duct. 
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Figure 34.- Phase velocities of mode n ■ 1 at various axial stations for 
q .667 at the inlet of a divergent duct. 
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Figure 35.- Phase velocities of mode n * -1 at various axial stations for 
q = 1.667 at the inlet of a divergent duct. 
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e n ■ 1 at various axial st. ions for 
et of a divergent dret. 


37 




